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ABSTRACT 

A  conparison  of  data  from  all  previous  energy  loss  experiments  per- 
formed at  the  Naval  Postgraduate  School  with  new  theoretical  predictions 
is  presented  in  this  thesis.  In  addition,  the  data  have  been  extended 
to  include  beryllium.  With  this  extension,  experiments  have  now  been 
conducted  on  materials  ranging  in  atomic  number  from  4  to  82.  Agreement 
between  experiment  and  theory  is  excellent  for  the  most  probable  energy 
loss.  However,  theoretical  values  for  the  half -widths  of  the  energy  loss 
distributions  generally  are  small  compared  to  experiment  for  thicknesses 
less  than  5  g/cm2  and  large  for  thicknesses  greater  than  5  g/cm2.  These 
experiments  were  conducted  in  the  energy  range  from  50  to  100  MeV,  with 

-i ,„i j_u,-  _i.„„-.. ,  .e— -_v,  r\  n    j-~  cr  r>  -,  /~~>2    mu„or,  ^.p-p^^-i-,^  r.,^. 

be  independent  of  atomic  number  or  incident  energy. 
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I.   INTRODUCTION 

Fast  electrons  which  traverse  a  layer  of  material  will  exit  with 
less  energy  than  they  possessed  on  entry.  The  energy  loss  distribution 
can  be  characterized  by  two  parameters,  the  most  probable  energy  loss 
(0  )  and  the  half -width  (HW)  ,  which  is  the  full  width  at  half  maximum 
of  the  distribution  curve.  These  quantities  are  the  basis  for  compari- 
son between  experimental  and  theoretical  results.  Losses  can  be 
attributed  to  two  major  effects,  ionization  and  radiation  (bremsstrah- 
lung) .  Ionization  losses  for  fast  electrons  have  been  calculated  by 
Landau  [1],  Williams  [2],  and  many  others.  Early  experiments  resulted 
in  good  agreement  with  theory  for  the  energy  loss,  but  these  results 
gave  a  half-width  which  was  greater  than  theory  at  energies  above  1  MrV. 
Blunck  and  Westphal  [3]  modified  previous  theoretical  calculations  by 
introducing  radiation  losses.  This  additional  energy  loss  is  super- 
imposed on  the  ionization  loss  and  results  in  a  larger  half -width  and 
enhancement  of  the  tail  of  the  distribution.  Other  losses  have  been 
discounted  as  unimportant  relative  to  ionization  and  radiation.  All  of 
the  theoretical  calculations  are  limited  to  "thin"  absorbers,  i.e., 
those  in  which  0  is  small  compared-  to  the  incident  electron  energy. 

Several  measurements  of  energy  loss  for  high  energy  electrons  through 
various  materials  have  been  conducted  at  the  Naval  Postgraduate  School. 
All  were  performed  within  an  energy  range  of  50  to  100  MeV.  These 
include  the  works  of  Bumiller,  Buskirk,  Dyer,  and  Miller  on  aluminum 
[4,5];  Goodwin  on  aluminum  and  copper  [6];  DeLeuil,  Raynis,  et  al  on 


aluminum,  copper,  and  lead  [7,8];  and  Mosbrooker  and  Sandquist  on  tin, 
lead,  and  gadolinium  [9].  With  the  extension  to  beryllium,  data  are 
now  available  on  various  materials  with  atomic  numbers  (Z)  between  4 
and  82. 

Experimental  consistency  with  past  experiments  was  insured  by 
rerunning  several  thicknesses  of  lead  absorbers.  Numerous  and  exten- 
sive checks  were  made  on  the  computer  program  (Appendix  D)  written  by 
Mosbrooker  and  Sandquist  [9]  to  insure  that  it  was  consistent  with 
Blunck  and  Westphal  [3]  theory  as  modified  by  the  beam  folding  tech- 
nique (Appendix  C)  which  accounts  for  the  energy  spread  of  the  incident 
beam. 

Comparison  of  all  experimental  data  taken  at  the  Naval  Postgraduate 
School  with  the  theoretical  results  of  Blunck  and  Westphal  [3]  is  pre- 
sented in  this  thesis.  The  data  on  beryllium  are  entirely  new;  the 
older  data  have  been  presented  before,  but  the  theoretical  values  given 
here  are  new  and  represent  calculations  made  with  the  computer  program 
of  Appendix  D. 


II.  THEORETICAL  CONSIDERATIONS 

The  Blunck  and  Westphal  theory  of  the  distribution  for  the  energy- 
loss  of  a  bean  of  monoenergetic  electrons  in  passing  through  a  layer 
of  absorbing  material  assumes  that  the  energy  loss  Q  is  small  compared 
to  the  initial  beam  energy,  E. .  Let  W(Q)dQ  be  the  probability  of  energy 
loss  between  Q  and  Q  +  dQ,  and  X  be  that  portion  of  the  loss  Q  due  to 
radiation.  Hence,  the  ionization  loss  is  Q  -  X.  Considering  these  two 
loss  processes,  the  probability  of  energy  loss  is 


Q 
W(Q)dQ  =  ;  WT(Q-X)W0(X)dXdQ 


(1) 


where  W  and  W_  are  the  energy  loss  distributions  for  ionization  and 
radiation  respectively. 

The  Landau  equation  [l],  as  modified  by  Blunck  and  Leisegang  [10] , 
is  used  for  the  energy  loss  distribution  due  to  ionization.  The  distri- 
bution is  expressed  as  a  function  of  Landau's  dimensionless  parameter 
A  as  follows: 


C  Y 

n'n 


WI(Q)dQ  =  $(A)  dX  -  Z 

n    b2  +  y  2 
n 


0—0 

where  X  =  ^—^    +  In 
aR 


1 
aR 


1.116 


exp 


(x  -  X  )2 
n 

dX 

(2) 

b2  +    y  2 
'n 

(3) 

The  terms  used  in  equations  (2)  and  (3)  are  defined  as  follows: 

C  ,  y  /  and  X  are  constants  given  by  Blunck  and  Leisegang  [10]  and 

are  used  to  fit  Landau's  distribution  to  a  sum  of  gaussian  functions. 
R  is  the  absorber  thickness  in  cm. 


The  quantity  "a"  is  a  function  of  the  atomic  number  Z,  the  atomic 
weight  A,  and  the  density  p  of  the  absorber;  and  3 (=v/c)  of  the 
electrons : 


=  0-154  Z  p 
A  32 


MeV/cm. 


The  quantity  b2  is  a  correction  to  Landau's  theory 
and  Leisegang  [10] : 


(4) 
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where  the  summation  is  over  the  ionization  potentials  of  the  atomic 
electrons,  and  N  is  the  number  of  electrons  with  ionization  potential 


V 


Q  is  the  average  energy  loss  due  to  ionization  (no  radiation)  for 


electrons  of  incident  energy  E. ,  and  is  given  by  Sternheimer  [11,12,13] 
as  follows: 


Q  = 


A  t 
s 


B  +  0.43  +  In  E. 


s2  -  C  -  ag(X1  -  log10p/mc) 


m 


MeV   (6) 


where  t  is  the  thickness  in  g/cm2,  and  the  constants  A  ,  B,  C,  X, ,  a  , 

and  m  are  parameters  of  the  absorber  material.  These  parameters  for 
tin,  lead,  copper,  beryllium,  aluminum,  and  various  other  materials  are 
listed  in  reference  [13] .  The  parameters  for  gadolinium  are  not  listed 
but  were  interpolated  by  Mosbrooker  and  Sanquist  [9] . 

For  W  ,  the  energy  loss  distribution  due  to  radiation,  Blunck  and 
Westphal  [3]  give: 


Ws(Q)dQ  - 


BoR(Q/E.)an  ^ 


(7) 
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where 


a  =  1.4  x  10 


-3  £Z 
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4/3  In  ^  +  1/9 
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cm 


-l 


(8) 


and  B  is  a  normalizing  factor  = 


r  (aR  +  1)   * 

The  distribution  of  total  energy  loss  according  to  Blunck  and 
Vfestphal  is  obtained  by  putting  equations  (7)  and  (2)  into  equation  (1) 
and  performing  the  required  integration.  The  result  is  the  energy  loss 
distribution  for  a  single  electron  of  incident  energy  E..  For  comparison 
of  theoretical  and  experimental  values,  this  distribution  function  was 
used,  with  corrections  to  account  for  the  finite  energy  width  of  the 
incident  electron  beam.  This  treatment  is  described  in  Section  IV  and 
Appendix  C. 
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Ill .   EXPERIMENTAL  PROCEDURE 

Beams  of  electrons  ranging  in  energy  from  53  to  95  MeV  were  obtained 
by  using  the  LINAC  of  the  Naval  Postgraduate  School.  Electrons  were 
elastically  scattered  from  the  accelerator  beam  through  90°  by  a  thin 
(1  mil)  aluminum  foil.  The  electrons  passed  through  the  varying 
absorber  thicknesses  and  were  energy  analyzed  by  a  120°  magnetic  spec- 
trometer described  by  Kenaston,  Luke  and  Sones  [14], 

The  general  experimental  .arrangement  wss  similar  to  that  used  by 
Mosbrooker  and  Sandquist  [9].   (See  Figure  1) .  The  energy  of  the  inci- 
dent electron  beam  was  determined  by  using  a  nuclear  magnetic  resonance 
probe  to  measure  the  field  of  the  deflection  magnet  of  the  LINAC.  The 
scattered  electrons  are  counted  on  a  ten  channel  plastic  scinfcillatoi 
counting  system  that  consists  of  ten  front  counters  and  a  single  backing 
counter  which  is  operated  in  coincidence  with  the  front  counters. 


SEM 

Al  Window 

Accelerator  Beam 


-Absorbers   ^Vacuum  Coupling 


Scattered  Beam 


Alternate  SEM  Location 


Figure  1.  Experimental  Arrangement 
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Beryllium  absorbers  of  thicknesses  ranging  from  0.738  g/cm2  to 
5.908  g/cm2  were  mounted  on  a  remotely  controlled  ladder  that  could  be 
moved  up  and  down  vertically  in  the  electron  beam.  This  allowed  us  to 
collect  data  en  one  absorber  thickness,  then  move  the  ladder  and  run  a 
different  thickness  without  having  to  enter  the  LINAC  end  station  or 
change  the  LIKAC  beam  in  any  way.  The  absorbers  were  placed  at  about 
3  cm  from  the  scattering  foil. 

A  lead  absorber  of  thickness  2.825  g/cm2  was  used  to  compare  our 
experimental  data  with  the  results  of  Mosbrooker  and  Sandquist  [9]  . 
The  results  were  in  agreement  within  experimental  accuracy. 

Since  the  electron  beam  scattered  into  the  absorbing  material  was 
not  monoenergetic ,  the  energy  distribution  was  measured  before  and 
after  passing  through  the  absorbers.  The  zero  thickness  absorber  dis- 
tributions found  by  passing  the  beam  only  through  the  scattering  foil 
had  half  widths  ranging  from  0.21  MeV  to  0.46  MeV,  or  an  energy  reso- 
lution of  about  0.45%.  Following  the  measurement  of  these  "zero  peaks", 
the  varying  absorber  thicknesses  were  moved  into  place  and  an  energy 
distribution  was  determined  for  each  absorber  thickness. 

Because  of  the  limited  quantity  of  beryllium  available  and  space 
requirements  in  the  target  chamber,  it  was  not  possible  to  place  all 
eight  absorber  thicknesses  on  one  ladder.  Therefore,  each  time  the 
beam  energy  was  changed  or  the  ladder  arrangement  was  changed,  a  new 
zero  peak  was  measured.  This  accounted  for  the  possibility  that  the 
character  of  the  beam  might  change. 

The  data  represent  the  number  of  scattered  electrons  detected  by 
the  coincidence  counting  system  at  the  exit  of  the  spectrometer  as  a 
function  of  the  energy  selected  by  the  spectrometer . 
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Normalization  was  controlled  by  integrating  the  current  frcm  a 
Secondary  Emission  Monitor.  This  instrument  was  mounted  both  up  and 
downstream  of  the  targets.  There  proved  to  be  no  significant  difference 
in  the  width  of  the  zero  peaks  due  to  the  position  of  the  monitor.  Each 
data  set  (ten  points)  represents  a  given  integrated  charge,  i.e.,  a 
certain  number  of  electrons  passing  through  the  scattering  foil. 
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IV.   TREATMENT  OF  DATA 

The  data  for  the  experiment  were  taken  with  the  ten  channel  counting 
system  as  described  in  Section  III.  To  reduce  the  data,  account  must  be 
taken  of  the  characteristics  of  the  counting  system.  These  include  the 
energy  spread  of  the  counters,  which  is  about  0.3%  between  counters  or 
3%  for  the  entire  systau.  Accordingly,  the  energy  seen  by  the  individual 
counters  must  be  determined  and  also  correction  must  be  made  for  the 
differing  efficiencies  of  the  scintillators  and  associated  electronics. 
A  computer  program  has  been  written  which  determines  the  correct  relative 
counter  energies  and  also  corrects  for  counter  efficiencies.  Background 
corrections,  which  were  quite  small,  were  determined  by  periodically 
removing  the  scatterer  and  absorbers  from  the  beam. 

The  reduced  data  were  then  plotted  and  from  these  plots  the  experi- 
mental most  probable  energy  losses  and  distribution  half -widths  were 
measured.  These  results  and  their  uncertainties  can  be  found  in  Table  I 
in  Appendix  A. 

Since  the  scattered  electron  beam  is  not  monoenergetic  and  the  Blunck 
and  Westphal  theory  is. based  on  a  monoenergetic  incident  beam,  to  compare 
theory  and  experiment  the  theoretical  model  must  be  modified.  This  modi- 
fication is  accomplished  by  representing  the  incident  distribution  as  a 
histogram,  each  bin  of  which  has  a  definite  energy,  and  then  unfolding 
this  distribution  into  the  theoretical  calculations  carried  out  by  an 
IBM  360/67  computer.  The  unfolding  technique  is  described  in  Appendix  C 
and  the  computer  program  used  is  in  Appendix  D. 
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The  computer  program  described  in  Appendix  D  also  calculates  the 
theoretical  half-widths  and  most  probable  energy  losses.  The  program 
normalizes  the  theoretical  curve  to  the  experimental  data  and  plots 
both  the  theoretical  curve  and  experimental  data.  These  plots  for 
beryllium  are  shown  in  Figures  [3]  through  [26]  in  Appendix  B.  The 
measurable  parameters  were  compared  and  the  final  results  are  tabulated 
in  Appendix  A. 

The  comparison  of  all  data  from  all  previous  energy  loss  experi- 
ments performed  at  the  Naval  Postgraduate  School  LINAC  was  accomplished 
by  normalizing  each  set  of  experimental  data  to  the  Blunck  and  Westphal 
theory  as  calculated  by  the  computer  program  described  in  Appendix  D. 
The  comparison  of  the  theory  parameters  so  obtained  to  the  experimental 
data  can  be  found  in  Tables  II  through  VI  in  Appendix  A. 
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V.   RESULTS  AND  CONCLUSIONS 

The  theoretical  predictions  and  experimental  results  for  the  most 
probable  energy  losses  and  half -widths  for  beryllium  can  be  found  in 
Table  I,  Appendix  A.  Our  theoretical  predictions  and  the  experimental 
results  from  previous  theses  completed  at  the  Naval  Postgraduate  School 
for  aluminum,  copper/  tin,  gadolinium  and  lead  can  be  found  in 
Tables  II  through  VT  respectively  in  Appendix  A.  The  theoretical  data 
shown  are  for  the  Blunck  and  Westphal  theory  as  modified  by  the  beam 
folding  technique.  In  all  tables  BW  refers  to  Blunck  and  Westphal  theory, 

0  represents  the  most  probable  energy  loss,  and  KW  is  the  distribution 
half -width.  The  column  headed  %  represents  the  percentage  by  which  theorry 
exceeds  experiment. 

Agreement  between  the  Blunck  and  Westphal  theory  and  our  experi- 
mental results  is  excellent  for  the  most  probable  energy  loss.  The 
average  difference  between  experiment  and  theory  is  less  than  0.5%. 
However/  agreement  for  half-width  is  good  only  for  thicknesses  less  than 

1  g/cm2  and  is  fair  for  thicknesses  between  4  to  6  g/cra2.  This  result 
in  half -width  comparison  is  in  contrast  to  previous  results  where  it 

was  concluded  that  the  Blunck  and  Westphal  theory  predicts  satisfactorily 
for  thin  targets  where  "thin"  was  defined  as  t(z)  1/3<13.5,  where  t  is 
the  absorber  thickness  (g/cm2)  and  z  is  the  absorber  atomic  number. 

The  observed  trend  for  HW  (experimental)  divided  by  HW  (theoretical) 
as  a  function  of  t  (absorber  thickness)  is  plotted  in  Figure  2. 
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FIGURE  2.  Trend  Plot  of  HW   /fiW..  vs  Thickness 

exp   th 


t(g/cm2) 


A  typical  comparison  between  experimental  and  theoretical  curves  for 
HW  plotted  as  a  function  of  t  (absorber  thickness)  is  presented  in  Figure  3, 


FIGURE  3.  Comparison  of  Experimental  and 

Theoretical  Half -widths  vs  Thickness 


t(g/cm2) 
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As  can  be  observed  from  these  plots  the  theoretical  half-width 
tends  to  be  smaller  than  experiment  until  one  gets  to  thick  absorbers 
(approximately  5  g/cm2)  at  which  point  it  gets  larger  than  experiment. 

We  have  considered  two  possible  effects  which  might  explain  why 
this  crossing  of  experiment  and  theory  takes  place.  One  effect  is  that 
this  trend  is  caused  by  an  assumption  that  is  made  in  the  theory.  This 
assumption  is  that  of  constant  incident  energy  (E.) ,  i.e.,  that  the 
energy  loss  is  small  compared  to  the  incident  energy  of  the  electrons. 
However,  in  thick  absorbers,  this  assumption  is  not  valid,  therefore 
the  theory  tends  to  overestimate  the  energy  loss  due  to  radiation 
(bremsstrahlung) .  This  overestimate  will  cause  an  increase  in  the  dis- 
tribution tail  and  hence  increase  the  distribution  half -width.  This 
trend  in  thick  absorbers  is  not  surprising  in  that  the  theory  was 
originally  designed  for  thin  absorbers.  The  fact  that  the  agreement  is 
excellent  for  0  while  not  for  half-width  is  also  not  surprising  since 
an  overestimate  of  this  type  would  not  affect  the  most  probable  energy 
loss  significantly.  Although  this  theoretical  inaccuracy  is  undoubtedly 
a  cause  of  some  discrepancy  in  agreement  with  experiment,  it  is  doubtful 
that  it  is  the  cause  of  all  the  disagreement. 

The  second  possibility  which  we  have  considered  affects  the  experi- 
mental results.  Due  to  multiple  scattering  in  thick  absorbers  the  possi- 
bility existed  that  the  size  of  the  spectrometer  opening  was  excluding 
the  low  energy  portion  of  the  distribution.  However,  calculations  of 
the  projected  mean  square  scattering  angle  as  outlined  in  Jackson  [15] 
showed  this  angle  to  be  approximately  0.30°  while  the  spectrometer 
opening  was  approximately  0.90°.  Therefore,  it  was  determined  that 
multiple  scattering  was  not  a  factor  in  the  disagreement  of  theory  and 
experiment. 
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Discrepancies  in  the  theoretical  predictions  from  our  program  and 
previous  results  of  Goodwin  [6]  and  DeLeuil  and  Raynis  [8]  were  noted. 
An  attempt  was  made  to  resolve  these  differences.  The  parameters  from 
Blunck  and  Westphal  [3]  were  inserted  into  our  computer  program  and 
without  the  unfolding  technique  described  in  Appendix  C  we  were  able  to 
duplicate  closely  the  Blunck  and  Westphal  distributions.  On  this  basis, 
it  is  assumed  that  our  theoretical  results  are  the  more  accurate. 

Contrary  to  previous  theses,  the  thickness  at  which  theory  and 
experiment  began  to  exhibit  a  noticeable  discrepancy  was  found  not  to 
be  dependent  on  the  atomic  number  of  the  material. 

It  should  be  mentioned  that  no  attempt  was  made  in  this  experiment 
to  observe  the  very  low  energy  tail  of  beryllium  (approximately  at  15 
MeV)  where  in  previous  experiments  secondary  effects  had  caused  a  peak 
in  the  absorber  distribution. 

In  light  of  what  has  been  found  in  this  thesis,  further  work  is 
justified  in  several  areas.  First  of  all,  the  low  energy  tails  of  the 
various  absorber  distributions  should  be  investigated  to  find  out  if 
secondary  effects  do  cause  a  low  energy  peak  in  some  distributions. 
Secondly ,  an  effort  should  be  made  to  improve  the  efficiency  of  the 
current  computer  program.  Finally,  an  attempt  should  be  made  to  extend 
the  theory  to  thick  absorbers  where  the  energy  loss  is  not  small  with 
respect  to  the  incident  energy,  and  also  to  include  the  effects  of 
secondary  processes. 
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APPENDIX  A  -  TABLES 


Table  I.  Energy  Loss  Distribution  Characteristics  of  Beryllium 


Be(z  = 

=  4) 

Qp(MeV) 

m  (MeV) 

E'.  (MeV) 

t(g/cm2) 

BW 

Experimental 

Q 
0 

BW 

Experimental 

O 

52.89 

0.742 

0.98 

0.98  ± 

.02 

-.30 

0.33 

0.36  ± 

.02 

-8.3 

1.479 

2.04 

2.00  ± 

.04 

+2.0 

0.53 

0.56  ± 

.04 

-5.3 

2.209 

3.05 

2.96  ± 

.04 

+3.0 

0.67 

0.77  ± 

.04 

-13.0 

2.961 

4.17 

4.09  ± 

.08 

+2.0 

0.90 

1.07  ± 

.07 

-15.9 

3.673 

5.29 

5.19  ± 

.05 

+2.0 

1.17 

1.25  ± 

.06 

-6.4 

4.415 

6.36 

6.24  ± 

.13 

+2.0 

1.40 

1.46  ± 

.09 

-4.1 

5.179 

7.46 

7.40  ± 

.08 

+  .80 

1.76 

1.76  ± 

.10 

0 

5.908 

8.69 

8.38  ± 

.12 

+3.6 

2.09 

2.19  ± 

.14 

-4.6 

74.78 

0.738 

1.02 

1.01  ± 

.03 

-1.0 

0.46 

0.47  ± 

.04 

-2.1 

1.479 

2.04 

2.01  ± 

.04 

+1.5 

0.56 

0.61  ± 

.04 

-8.2 

o    inn 

Z  ,/L\JZJ 

->     n 

■a    r\A    + 

~>  •  U'l      X 

.  \j~> 

+2.3 

C\     -7/1 

n    TQ    x 

O    .     1    ^         — 

.05 

-6.3 

2.941 

4.15 

4.15  ± 

.14 

0 

0.96 

1.15  ± 

.U8 

-16.5 

3.673 

5.29 

5.15  ± 

.06 

+2.7 

1.17 

1.39   ± 

.06 

-15.8 

4.435 

6.39 

6.28  ± 

.10 

+1.7 

1.46 

1.61  ± 

.10 

-9.3 

5.179 

7.46 

7.42  ± 

.13 

+  .50 

1.76 

1.83   ± 

.08 

-3.5 

5.908 

8.68 

8.51  ± 

.09 

+2.0 

2.12 

2.24   ± 

.11 

-5.3 

94.64 

0.738 

1.02 

1.05  ± 

.04 

-2.7 

0.58 

fi«37  ± 

.05 

-1.7 

1.479 

2.04 

2.06  ± 

.04 

-1.0 

0.68 

0.66  ± 

.04 

+3.0 

2.209 

3.11 

3.09   ± 

.04 

+  .60 

0.83 

0.84  ± 

.06 

-1.2 

2.941 

4.15 

4.06  ± 

.06 

+2.7 

1.06 

1.22  ± 

.08 

-13.1 

3.673 

5.29 

5.19  ± 

.06 

+1.9 

1.23 

1.35  ± 

.08 

-8.9 

4.435 

6.39 

6.30  ± 

.12 

+1.4 

1.53 

1.58  ± 

.12 

-3.2 

5.179 

7.61 

7.47  ± 

.12 

+1.9 

1.82 

1.93  ± 

.10 

-5.7 

5.908 

8.69 

8.63  ± 

.14 

+  .60 

2.16 

2.29  ± 

.11 

-5.7 
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Table  II.  Energy  Loss  Distribution  Characteristics  of  Aluminum 


A1(Z  = 

=  14) 

Q  (MeV) 

HW(MeV) 

Ei  (MeV) 

t(g/cm2) 

BW 

Experimental 

g. 
"5 

BW 

Experimental 

Q. 

53.57 

0.730 
1.441 
2.146 
2.859 

1.05 
2.12 
3.22 
4.29 

1.03 
.2.16 
3.18 
4.21 

+1.9 
-1.9 
+1.2 
+1.9 

0.44 
0.62 
0.88 
1.20 

0.49 
0.75 
1.02 
1.32 

-10.0 

-17.  Oj 

-14.01 

-9.1 

53.44 

5.574 

8.86 

9.24  ±  .13 

-4.1 

3.12 

4.60  ±  .30 

-32.0 

74.63 

0.730 
1.441 
2.146 
2.859 

1.05 
2.12 
3.22 

4.37 

1.04 
2.18 
3.17 
4.26 

+1.0 
-2.7 
+1.6 
+2.6 

0.54 
0.72 
0.95 
1.25 

0.51 
0.91 
1.08 
1.38 

+5.9 
-21.0 
-12.0 

-9.5 

74.47 

5.574 

8.86 

9.33  +  .25 

-5.0 

3.07 

5.7  ±  1.2 

-45.0 

96.93 

0.730 
1.441 
2.146 
2.859 

1.10 
2.16 
3.22 
4.37 

1.17 
2.16 
3  24 
4.32 

-6.0 

-0.6 
+1.2 

0.64 
0.81 
1.03 
1.31 

0.66 
1.01 
1.19 
1.48 

-3.0 
-20.0 
-13.0 
-11.0 
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Table  III.  Energy  Loss  Distribution  Characteristics  of  Copper 


cu(z  = 

=  29) 

C)  (MeV) 

HW(MeV) 

E^  (MeV) 

t(g/cm2) 

BW 

Experimental 

o. 

BW 

Experimental 

Q. 

52.84 

0.711 
1.423 
2.134 
2.845 

0.92 
1.89 
2.90 
3.87 

0.94 
1.87 
2.87 
3.80 

-2.1 
+1.1 
+1.0 
+1.8 

0.46 
0.69 
1.04 
1.54 

0.52 
0.84 
1.16 
1.75 

-12.0 

-i8.o; 

-10.  Oi 
-12.  Oj 

53.63 

4.295 

6.10 

6.23  ±  .13 

-2.1 

3.70 

4.25  ±  1.20 

1 

i 

-13.0 

74.76 

0.711 
1.423 
2.134 
2.845 

0.92 
1.89 
2.90 
3.96 

0.92 
1.91 
2.92 
3.98 

-1.0 
-0.7 
-0.5 

0.56 
0.78 
1.12 
1.61 

0.66 
1.07 
1.43 
1.98 

-15.0 
-27.0 
-22.0 
-19.0 

76.02 

4.295 
5.726 

6.10 
8.65 

8.38  ±  .25 

-2.9 
+3.2 

O   TO 
_)  •  /  O 

10.37 

4.55  ±  .50 
10.1  ±  2.0 

-17.0 
+3.6 

94.30 

0.711 
1.423 
2.134 
2.845 

0.91 
1.94 
2.96 
3.99 

0.95 
1.94 
2.93 
4.00 

-4.2 

+1.0 
-0.2 

0.67 
0.87 
1.22 
1.71 

0.74 
0.99 
1.47 
1.75 

-9.5 
-12.0 
-17.  C 

-2.3 

96.66 

5.726 

8.82 

9.20  ±  .50 

-4.1 

10.50 

16.4  ±  4.2 

-37.0 
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Table  IV.  Energy  Loss  Distribution  Characteristics  of  Tin 


Sn(Z  = 

=  50) 

Q  (MeV) 

HW(MeV) 

Ej^MeV) 

t(g/cm2) 

BW 

Experimental 

ft, 
t5 

BW 

Experimental 

a 
o 

52.53 

1.485 
2.970 
4.455 

1.83 
3.84 
6.56 

1.78  ±  .05 
3.86  ±  .15 
6.35  ±  .30 

+2,8 
-0.5 
+3.3 

0.75 

3.03 

16.70 

0.78  ±  .10 

3.02  ±  .40 

11.30  ±  2.0 

-3.8 

+0.3 

+48.0 

75.00 

1.485 
2.970 

4.455 

1.81 
3.93 
6.56 

1.77  ±  .07 
3.86  ±  .15 
6.65  ±  .30 

+9.0 
+1.8 
-1.4 

0.73 

2.56 

16.69 

0.83  ±  .10 

2.96  ±  .42 

15.50  ±  2.0 

-12.0 

-14.0 

+7.7 

94.40 

1.485 
2.970 
4.455 

1.88 
3.93 
6.56 

1.87  ±  .10 
3.96  ±  .20 
7.00  ±  1.0 

+0.5 
-0.8 
-6.3 

0.83 

2.63 

16.75 

0.88  ±  .15 

2.72'  ±  .45 

16.20  ±  4.0 

-5.7 
-3.3 
+3.4 

Table  V.  Energy  Loss  Distribution  Characteristics  of  Gadolinrum 


Gd(z  ■- 

=  64) 

Q^MeV) 

HW(MeV) 

EL|_  (MsV) 

t(g/cm2) 

BW 

Experimental 

Q. 

BW 

Experimental 

% 

53.18 

0.814 

0.97 

0.92  ±  .08 

+5.4 

0.46 

0.52  ± 

.07 

-12.0 

52.80 

1.610 
3.221 

1.94 
4.24 

1.92  ±  .10 
4.16  ±  .30 

+1.0 
+1.9 

0.91 
5.52 

1.00  ± 
5.75  ± 

.15 
.72 

-8.9 
-4.0 

75.00 

0.814 
1.610 
3.221 

0.97 
1.97 
4.24 

0.97  ±  .05 
1.92  ±  .10 
4.09  ±  .30 

+2.6 
+3.7 

0.51 
0.96 
5.55 

0.59  ± 
1.00  ± 
5.80  ± 

.08 
.15 
.80 

-13.0 

-4.0 
-4.2 

91.80 

0.814 
1.610 
3.221 

0.97 
1.97 
4.20 

0.97  ±  .05 
1.90  ±  .10 
4.13  ±  .40 

+3.7 
+1.7 

0.60 
0.96 
5.56 

0.64  ± 
1.10  ± 
8.10  ± 

.08 
.18 
.95 

-6.2 
-13.0 
-31.0 
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Table  VI.  Energy  Loss  Distribution  Characteristics  of  Lead 


Pb(Z 

=  82) 

Q  (MeV) 

HW(MeV) 

Ej_  (MeV) 

t(g/cm2) 

BW 

Experimental 

BW 

Experimental 

Q. 

53.85 

0.706 

0.79 

0.81  ± 

.05 

-2.5 

0.54 

0.60  ± 

.06 

-10.0 

1.412 

1.66 

1.45  ± 

.08 

+14.0 

0.97 

1.23  ± 

.12 

-21.0 

2.118 

2.55 

2.64  ± 

.10 

-3.4 

2.11 

2.65  ± 

.25 

-20.0 

2.825 

3.57 

3.66  ± 

.12* 

-2.5 

6.45 

6.25  ± 

.50* 

+3.2 

74.74 

0.706 

0.83 

0.82  ± 

.05 

+1.2 

0.78 

0.82  ± 

.05 

-4.9 

1.412 

1.68 

1.80  ± 

.10 

-6.7 

1.23 

1.47  ± 

.07 

-16.0 

2.118 

2.66 

2.59  ± 

.12 

+2.7 

2.37 

2.24  ± 

.19 

+5.8 

2.825 

3.66 

3.70  ± 

.15* 

-1.1 

6.52 

6.25  ± 

.85* 

+4.3 

91.37 

0.706 

0.90 

0.87  ± 

.07 

+3.5 

0.93 

0.78  ± 

.08 

+19.0 

1.412 

1.80 

1.77  ± 

.10 

+1.7 

1.37 

1.52  ± 

.19 

-9.8 

2.118 

2.70 

2.75  ± 

.15 

-1.8 

2.53 

2.74  ± 

.64 

-7.7 

.  - 

2.825 

3.69 

3.84  ± 

.22 

-3.9 

6.95 

8.80  ± 

.70 

-22.0 

*Experiment  data  average  of  DeLeuil/Raynis  and  Mosbrooker/Sandquist 
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APPENDIX  C  -  BEAM  FOLDING  TECHNIQUE 

The  theory  of  Blunck  and  Vfestphal  [3]  assumes  that  the  beam  of 
electrons  striking  the  absorber  is  monoenergetic.  The  beam  of  elec- 
trons produced  by  the  Naval  Postgraduate  School  LINAC,  or  any  other  linear 
accelerator,  for  that  matter,  has  a  finite  energy  spread  about  a  maximum 
or  most  probable  energy  point.  The  fact  that  monoenergetic  electrons 
are  not  available  to  strike  the  absorber  must  be  taken  into  account  in 
computing  theoretical  predictions  if  a  meaningful  comparison  with 
experimental  results  is  to  be  made.  This  has  been  done  here  by  a 
technique  termed  "beam  folding". 

Beam  folding  is  accomplished  by  a  number  of  well  defined  steps. 
Energy  distribution  carves  may  be  approximated  by  histograms.  These 
histograms  consist  of  a  series  of  "bins"  of  area  W(Q)AQ,  where  Q  is 
the  energy  loss  and  W(Q)AQ  is  the  probability  of  loss  between  Q  and 
Q  +  AQ.  Thus,  each  bin  has  an  "address",  Q,  on  an  energy  coordinate 
scale.  To  accomplish  beam  folding,  the  beam  distribution  must  be  known. 
This  is  observed  experimentally  and  approximated  in  the  computer  by  a 
gaussian  curve  of  appropriate  half-^idth.  The  energy  at  which  the 
maximum  occurs  is  established  by  the  energy  of  the  beam  and  the  magni- 
tude of  energy  loss  incurred  by  elastic  scattering  of  the  beam  as  it 
impinges  on  the  thin  aluminum  scattering  foil  prior  to  striking  the 
absorber.  The  width  AQ  is  then  selected  for  the  predicted  distribution. 
This  must  be  a  small,  but  finite  number  where  numerical  techniques  are 
to  be  used.  This  same  width  is  used  to  break  the  beam  distribution  into 
a  histogram.  This  is  illustrated  in  Figure  30  (a)  .  The  reason  the  sa~2 
width  is  used  for  both  distributions  is  to  facilitate  computer  prograimiing . 
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Each  bin  of  the  beam  distribution  is  now  treated  as  a  monoenergetic 
beam  with  energy  commensurate  with  its  center  location  on  the  energy 
scale.  Each  bin  also  has  a  definite  magnitude,  or  weight,  with  the  magni- 
tude of  the  center  bin  being  unity.  The  formulae  of  Blunck  and  Vfestphal 
is  now  applied  to  each  of  these  beams  and  an  absorber  distribution  curve 
results  for  each,  with  a  maximum  amplitude  proportional  to  the  height  of 
the  appropriate  beam  distribution  histogram.  Each  of  the  absorber  dis- 
tribution curves  thus  obtained  may  be  thought  of  as  being  plotted  and 
added  to  previously  determined  curves,  using  the  energy  of  the  electrons 
as  a  correlation  point,  as  depicted  in  Figures  30(b) ,  (c)  and  (d) .  This 
is  accomplished  on  the  computer  by  adding  the  contents  of  each  bin  of  the 
same  address  and  plotting  the  cumulative  total.  The  bin  of  address  E  is 
depicted  in  Figures  30(b),  (c)  and  (d) .  Note  that  the  beam  distribution 
for  histogram  E4  is  centered  over  the  maximum  point  of  the  distribution. 
This  is  important  as  a  false  picture  could  easily  be  presented  if  the 
histogram  were  not  symmetrical  as  the  distribution  would  then  appear 
skewed. 
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APPENDIX  D  -  CCMPUTER  PROGRAM 


C  PURPOSE:  THIS  PROGRAM  ACCEPTS  PARAMETERS  FROM  ENERGY  LOSS 

C    EXPERIMENTS  CONDUCTED  ON  THE  LlNAC  AND  WILL  COMPUTE 

C    THEORETICAL  HALFwIOTh  AND  MOST  PROBABLE  ENERGY  LuSS  AND 

C    WI LL  PLJi  BOTH  THE  EXPERIMENTAL  AND  THEORETICAL  ENERGY 

C    DISTRIBUTION  FOR  THE  EXPERIMENT, 

C 

C  SCOPE:  THE  PROGRAM  IS  SET  UP  TO  MAKE  COMPUTATIONS  RELATIVE 

C    TO  EXPERIMENTATION  WITH  ALUMI NUM. COPPER, T I N, GADOL IN IUM 

C    LEM  AND  BERYLLIUM. 

C 

C  METHOD:  THE  PROGRAM  PLOTS  A  THEORY  CURVE  BASED  ON  THE  FIT 

C    OF  GAUSSIAN  CURVES  PERFORMED  BY  BLUNCK  AND  WES1PHAI.  FOR 

C    THf  PREDICTED  ENERGY  DISTRIBUTION  CORVES,   THE  INTEGRA- 

C    TION  OVE«  ENERGY  LOSS  INDICATED  IN  The  AFOREMENTIONED 

C    PAPER  IS  PERFORMED  ON  THE  COMPUTER  USING  3?  POIN1  GAUSS 

C    QUADRATURE.   TWO  GRAPHS  *RE  PLOTTED;  ONE  REPRESENTS  1  HE 

C    B-VJ  PREDICTION  WITH  A  ''.ONENERGET  IC  BEAM  OF  ELECTRONS  AND 

C    TE  !-  OTHER,  WHICH  ALSC  HAS  EXPERIMENTAL  RESULTS  THEREON, 

C    REPRESENTS  TH^  SAME  THEORY  WITH  A  POLYENbRGET IC  BEAM 

C    DISTRIBUTION  FOLDED  IN. 

C 

C  THIS  PROGRAM  IS  WRITTEN  USING  FORTRAN  IV  LANGUAGE  AND  THE 

C    NPS  COMPUTER  FACILITY  PLOTTING  PACKAGE.   TOTAL  RUNNING 

C    TIME  IS  A  FUNCTION  OF  THICKNESS,  ATOMIC  NUMBERi  AND 

C    HALFVIUTH. 

C 

I  vfLlCl  T  REAL* 8 ( A-H.O-Z ) 
C 

C  ESTABLISH  TITLFS  FOR  PLOTS. 
C 

REAL*8  Tl<4) 

CATA  Tl/«        ENERGY  DISTRIBUTION       •  / 

REAL*8  T3(4) 

CATA    T.i/»  ALUMINUM    ABSORBER  •  / 

REAL* 8     TMM 

RATA    T4/'T:  GM/CM  LI  :  MEV  / 

REAL*8     T5(4) 

DATA  T5/1         UNFOLDED   THEORY         «/ 

REAL*8  T6(4) 

CATA    Tft/'THEORY:  EXPERIMENTAL    DATA:  XV 

REAL-  8     77(4) 

CATA  7  7/'  LEAD  ABSORBER  •  / 

REAL* 8  T8(4) 

i:\TA  T3/«         TIN  ABSORBER  '  / 

R  t  \  L  *  8  ■  T  9  (  4  ) 

RATA  79/ «         GADOLINIUM  ABSORBER        •/ 

REAL*8  71G<4) 

RATA    T1U/'  COPPER  -ABSORBER  V 

REAL* 8    7  1KM 

CATA  Til/1        BERYLLIUM  ABSORBER        « / 

DIMENSION  E(50) ,WF( 5  ) »WQT(3  I  ),HW( 30  0 , XXC30  )) » 
iWUF(3cl  )  ,XUF(3G?M  ,EXX(?A'0<)  ,EXY{2v;G)  ,  IXY(2DG),72(4) 

EXTERN AL  F C T , BEAM, W I  ON , WRAD 

DC12C  U  J  =  l,5 

E  (  J  )  =  0 
12  00  W F ( J ) = 

DC12v-l  J  =1,200 

EXX( J)=( 

EXY(J)=0 
12CI  IXY(J)=, 

D012D2  J  =  l,3i  C 

WGT ( J)=0 

;  I  •  (  J  )  - 

X  X  (  J  )  = 

W  U  i-  (  J  )  =  ) 
1202  XUF( J)=  ) 
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c 
c 
c 
c 


c 
c 

c 
c 
c 

c 
c 
c 


READ  IN  PARAMETERS  FOR  EXPERIMENT, 
EB       BEAM  ENERGY 


D 


21    FORMAT (/,2X,'     EB:     • tF7.2»'     Z:     '^4.0,'     A:     •  ,F8.3f 
}'     T:     '  ,F7.4,»     .JZtl'u:     SF5.3,1     CTSNOR:     • ,F6.1, 
2'    HT :     '  , Fb.3»/  ) 

IVlOfc    ZERO    PEAK    INTO    "BINS11    OF     WIDTH    DELX.       BIN    WIDTH 
I N C R E A S F S    W I T H    T A R OFT    THICKNESS. 

DELX  =  .i    V-T 

IFtT.Ll . 1, )    DELX=.06*T 

QAVE  IS  THI  RECOIL  LOSS  FF  ELECTRONS  STRIKING  ALUMINUM 
NUCLEI  1  OF  FHE  SCATTERING  FOIL.   IONIZATION  LOSS  IN 
SCATTEREP  IS  NEGLECTED. 

QAVE=EB*EB/< 931. 4*27. ) 

CORR  IS'THf  CORRECTION  TO  RE  ADDED  TO  EXPERIMENTAL  ENERGY 
TO  CORRECT  FOR  SPECTROMETER  ERROR. 

CORR=QZERO-QAVE 
12  FORMAT (5(1X,F8.3,1X,I6) ) 

READ  IN  UP  TO  2  £  EXPERIMENTAL  POINTS;  EXX  IS  ENERGY  (MEV) 
4ND  IXY  IS  COUNTS.   THERE  MUST  RE  4(J  DATA  CARDS  ARE  THE 

LAST  FXX  MUST  BE  ZERO. 

RE AD ( 5 , 1 2 )  ( E X  X ( J  )  ,  I  X  Y ( J )  , J  - 1  , 200  ) 

ESTABLISH  PARAMETERS  FOR  SCATTERING  FOIL. 

FOIL-.(  }686 
ZF0IL=1 3. 
KSIG=0 
DO  1  J -1,2.; 
I  EXY( J)=IXY( J) 

EI  IS  VALUE  OF  BEAM  ENERGY  INCIDENT  ON  TARGET. 

EI=EB-QAVE 
TE  4P=E  l 

THE  FOLLOWING  PORTION  OF  THE  PROGRAM  UNFOLDS  BEAM  DISTRI- 
BUTION.   ICENT  IS  AN  INTEGER  VALUE  ASSIGNED  TF  THE 
HISTOGRAM  CENTERED  OVER  THE  PEAK  LIE  THt  BEAM  DISTRIBU- 
TION.  EACH  HISTOGRAM  IS  DENOTED  BY  ITS  ENERGY  F  AND 
MAGNITUDE  RELATIVE  TO  THE  PEAK  OF  THE  DISTRIBUTION,  WF. 

ICENT=(HT+DELX ) /DFLX 
X=EI-ICENT*DELX 
202  X=X+DELX 

IF(X.Gr. (EI+HT  )  )  GOTO  1003 
M=M+1 
E  ( M )  =  X 

WF(M)=BEAM(EI,HT,X) 
GOT 02  ! 
1   13  WRIT  (  >  ,  :\?  )  DAVE  tHTT  FOIL,  ZF  OIL  t  EI 
1   •   L=l, 
DX=j 
N=  ) 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


1 


USING  EACH  HISTOGRAM  AS  A  "BEAM",  THE  FOLLOWING  PORTION  OF 
THE  PROGRAM  COMPUTES  AN  ENERGY  0  1  S  1  R  I  BUT  1  ON  FOR  EACH 
OF  THESE  "BEAMS"'  AND  STORES  THEM  IN  "BINS"  LABELED  BY 
ENERGY,  XX.   THE  CUMULATIVE  AMOUNT  OF  ELECTRON  COUNTS  IN 
THE  JIN  IS  WOT.   Q  IS  THE  ENERGY  LOSS  OVER  WHICH  THE  B-W 
PREDICTION  IS  INTEGRATED  AND  DX  IS  THE  AMUUNT  BY  WHICH 
Q  CHANGES  fO:-     EACH  SUCCESSIVE  INTEGRATION,   INTEGRATION 
IS  CARRIED  OUT  IN  SUBROUTINE  DQG32. 


X=DX+DELX 

F(Z.LT, 10. )  GO  TO  333 

=  3.0*T*(Z/13.  )**.  S.  -nx 


3  33 
334 


li  04 


u 

I 

Q 

GG  TO  334 


C 

c 
c 
c 
c 

r 


=3.0*T*<Z/4. ) **.5  -OX 
F(Q.LT.EI)  GO  TO  L0u4 


-CTOlf   5 


MQ. 
I  =  t( 
L  =  .  ■ 
ALIM 
ALL 
Ml. 
ALL 
G=ZA 
Q  T  ( L 
X(  LJ 
=  L  +  0 


LI 
L) 


U ) GOTOl 006 

j:  i*q 

R*0 


I)  3G 
DOG 
DQG 
+  ZB 
+  N- 
N-l 
-1 


.rCT, 

,FCT, 


32U., XL 

32 ( XL, VALI 

32  (  VAL  I  -V) 

fZC 

1  )  =  WQT  (  L  +N-1  )  +  «'Q 

)=E(L)-Q 


, FCT,ZA, LI 


Z  3  , 
ZC, 


El 
El 


,  Q ,  T  ,  A ,  Z  ) 

t  Q  ♦  T  t  A ,  Z ) 


t  i  , 


A ,  Z  ) 


WF  (  L) 


THE  DISTI  IPUT  UN  FOR 
Of  MTfc  R  OF  THE  Br 
WUF  \ND  *UE  AS  THE 
RESPECT  IV ELY. 


THE  HISTOGRAM  CORDESPONDI NG  TO  THE 
DISTRIBUTION  IS  P., MOVED  AND  STORED 
UNFOLDED  COUNTS  AOD  ENERGY 


I    ! 


C 

c 
c 
c 
c 


I F ( L .NF. ICE  NT) GOTO  10 J5 
hUF(N)  --.;  j 
xijf  (N>  =EI-Q 
K=N 

GOTO    1  ■.   o 

CCNT  INUE 

DO    1717     1=1,7 


ELIMINATE  LGK  NUMBERED  "BIOS"  ".S 
LEO  \ND  ARE  NOT  REPRESENTATIVE 
DISTRILUJT  ICO. 


THFY  ABE  ONLY  PARTLY  F I L- 
OF  THE  CUMULATIVE 


C 

c 
c 
c 
c 


feUF(  I  )  =  njF( 

17  17  WQT(.I)=WQT( 

DO  6  1=1,4 

6  T  U   I  )  =  T  6  (  I  ) 

GOTO  779 

77  7  DG779  N=1,K 


3) 
3) 


0 

c 
c 


PROCESS  THE  UNFOLDED  THEORY  CURVE 
TO;  ZERO  ALLOWS  PLOTTING  OF  THIS 
EXPERIMENTAL  DATA. 


w'U  (0) =WUF( N) 

XX(N)=XUF(N) 

KSIG=1 

J=K 

CONTINUE 
EXX(1)=0 
DO  5    1=1,4 
T  3  (  I  )  =  T  5  (  I  ) 

rMx=. 

0  10  J  7 


SETTING  EXX(l)  EQUAL 
THEORY  CURVE  W  ITHOUT  A  l\ 


7  76 


5 
779 


1=3, J 


DETERMINE  THE  .'-  A  X I  MUM  PL 
AND  THUS  ESTABLISH  THE 


INI   F  THE  0ISTRI3UT ION, 
MOST  PROBABLE  ENERGY  LOSS,   .  < 
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IF(  WQT(N).LE.WQT(N-1  )  .AND.WQ1  (N-l  )  .  GB  .  WOT  (  N-2  )  ) 
1GOT0100  8 

GOTO  LC07 
LC08  IF(WQT(N-l).LT.WGTMX)GOTO  L007 

WQTMX=WQT ( N-L) 

E I  =  T  E  M  P 

3P=E I -XX (N-l ) 
1. 07  CONTINUE 

NORMALIZE  ALL  POINTS  OF  DISTRIBUTION  AND  STORE  IN  HW. 


C 

C 

c 


DO  1009  M=1,J 
HW(M)=WQT(M)/WQTMX 
1CC9  CONTINUE 
J=J-1 
HWI=G 

DETERMINE  HALFWIDTH  OF  THE  THEORETICAL  DISTRIBUTION!  HWA. 

D  G  11   H  =  I ,  J 

IF(HW< M  +  l) .GT. ,5.AND,HW(M),LT..5)HWl=XX(M)+(.5- 
1HW(M) ) *DELX/{HW( M+l )-HW(M) } 
1 F ( M. EQ ♦ J • A  NO. HW 1  * E Q.G ) GGTOl 1 9 
GOT0121 
119  WRITE (6, ?  ) 
HWA  =  U 
GOTO] 01 0 
20  FORM/\T(/,2X,»NO  HALFWIDTH  OBTAINED') 
121  IF(HW( M) .G1 . .5.AND.HW( M+l ) . LT. . 5 ) HWA=XX ( M ) +  (HW< M) 
1-.  5  )>>  DE  L  X/  (  HW  (  M  )  -H  H  (  M  +  l  )  )  -HW  I 
K  10  CONTINUE 

ESTABLISH  CUTOFF  SIGNAL  FOR  PLOTTING  SUBROUTINE  AND 
ESTABLISH  CORRECT  GRAPH  TITLES  FOR  ADSORBER  USCD. 

V/f  J  +  1  )  -  | 

IF( Z.GT  .4.1  ) GUT  01 3 
DO  2  1  =  1 , 4 

2  r  2  (  I  )  =  T  1  1  (  I  ) 
GOTO101 1 

13  IF(Z.GT.L3.1)GOT029 
DO  3  1=1,4 

3  T  2  {  I  )  =  T  3  (  I  ) 
GOTO! 01  1 

2^  IF(Z.GT.29.1)G0T05G 

DO  101=1,4 
10  T2( I)=T10( I ) 

GOT 01 01  1 

IF(Z.GT.50.1)G0T064 

DO  8  1=1,4 
8  T2( I)=T8  ( I ) 

GO TO 10 11 
64  IF(Z.GT.64.1)G0T082 


DC 


1  =  1,4 


9  T  2  (  I  1  =  T  9  (  I  ) 
GO  TO  101 1 
82  DO  7  1=1,4 
7  T  2  {  I )  =  T  7  (I) 

SEND  DATA  TO  PLOTTING  SUBROUTINE  FOR  PROCESSING  AND 
PRINT  RESULTS  OF  COMPUTATIONS  IN  TABULAR  FORM, 

1.11  CALL  GRAPH ( XX,WQT,EXX,EXY,T4 ,T3 ,T2 , Tl ,HWA, EI , T, 
1CORR,CTSNOR,OP ) 
WRITE (6,22)  OP, HO  A, T,Z,EI 
IF(KSIG.EQ.O)  GOTO  777 
1 1  FORMAT ( 7 Fl  .  ) 


?2    F   M A  T ( /  ,  2  X  ,  « 

1  '  T(GM/S 3-C  I  : 

2  ,  F  6  .  2  ,  /  ) 

SI 

END 


'  ,  F  7 . 

F6.4, • 


'•+  » 


Z: 


HAl  F  W  IOTH:  •  , F6.3, 

•  ,F4.l  , •  INC!  DENT  f    NY:  ' 
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PURPO 

IN 

I 

A 
B 
C 
Y 
1  ( 

c 

Y 

1  ( 
C 
Y 

1  ( 
C 
Y 

.1  ( 
C 
Y 

1  ( 
C 
Y 

1  ( 
C 
Y 

1  ( 
C 
Y 

1  ( 
r 


SUBROUTINE  DQG32 (XL ,XU t FCT , Y, E  I  ,Q tT, D, Z ) 

ERFORMS  THE  INTEGRATION  REQUIRED 


S3:  THIS  SUBROUTINE  P 
THE  B-W  THEORY. 


1  ( 

C 
Y 

1  ( 
C 
Y 

1  ( 

C 

v 

1  ( 

C 
Y 

1  ( 
C 
Y 

1  ( 
C 
Y 

1  ( 
C 
Y 

1  < 
R 
E 


KPLICI 

=  .5D  >* 
=XU-XL 
=  .49  Jo 
=. 350  9 
FCT  (  (  A 
=.492  8 
=Y+.81 
FCT(  (  * 
=.4623 
=  Y  +  . 1  I 
FCT(  (/' 
■=.4  674 
=  Y+.  17 
FCT( (A 
=.4431 
=Y+.21 
FCT( (A 
=  .  4  2  4  t 
= Y  +  .  2  '  5 
FCT( (A 
=.397? 
=Y+. 29 
FCTU  . 
=  .  366'. 
=  Y  +  .  3  2 
FCT( (A 
=  ,  ^  !  r> 
=Y+.36 
F  CT ( {  * 
=.2938 
=Y+. 3  9 
F  r  t  (  (  A 
=  .253<- 
=Y+.41 
FCT( (A 
,=  .2lOt 
'=Y  +  .  ■'*  ? 
F  C  T  (  (  A 
=.1659 
=  Y  +  .  4  5 
FCT{ (A 
.=  .  119  6 
=  Y+.46.' 
FCT  (  (A 
=.7223 
=Y+.47 
FCT(  (  * 
= .  2  4 1  S 
=B*(Y+ 
FCT( ( ft 
ETURN 
NO 


T    REAL*8( A-H,n-Z) 
(XU+XL) 


v- 


3  19  3  ■? 
3  5  C47 
4-C)  ,/,i) 
0  5  7 '3 5 7 7 
3  719  7  3  6 
4-C)  ,Z,D 
8  112  779 
696  .*326 
4-C)  f  ZtO 
5303 796 
136  9  314 
*  C )  ,  Z ,  0 
6 3  5  7  7 : 8 
4 1 794 9u 
4-C)  ,Z,D 
8  38  t  '■'• 
499U2  )6 
+C) , 7 ,0 
4189 798 
3  -'♦  7  ■  ■  m  6  7 
4-  C )  ,  Z  ,  0 
91J59 j 7 
9111113 

+  C )  ,  z  t  o 

?  2  I :     v  > 

17  3  - 


578  7 

5  969478 
4-C)  ,Z,0 
499  5  4 A  6 

6  559 '  21 
4-C)  rZ  ,D 

7  5  6  3 i  '6 
82  6  46  5 
4-C)  ,  7,D 
343  114 
586939 3 
4-C)  ,Zt  0 
\  36  !  1 1 2 
9221995 
4-C)  ,7,0 

3 -J  8  791 
8193  6UU 
4-C)  ,  7  ,  0 
3832843 
.  4  3  2  7  .   5 

4"  C  )  ,  /  ,  L) 


47  4  U 
3  5  J  4 
t  T ,  E 
2634 

->  **■'■>  2 
,T,E 
3  7  3  3 
34  6  3 
,T,E 
3369 
i"i  a  t_»  1 

i  T  ,  E 
30  2  6 
1111 
,T,E 
33  84 
3113 
,  T  ,  E 

39  71 
3  9  2  6 
,T,E 

I  4  4 
88  18 
f  T,E 
5 1  7 
>  4  4  2 
f  T,E 
,  3  8  1 
)  3  3  3 
i  T,t 
6114 
1  3  47 
,T,E 
5317 

220 
,T,E 
L.  .63 
47  8  8 
,T,E 
6  :op, 

40  4 'J 
,  T,E 
3  9  82 
3  9  6  3 
,  T  ,  E 

3  6  9 1 

4  4  2  5 
J,b 


7  800 

8  30- 

1  ,  0  ) 
1700 
835D 
1,0) 

2  2  DC 

l  e  3  < 

I  ,  Q ) 
840v 
.717 
1,0) 

.  6o: 

3  34( 
I,Q) 
99  i 
34  3  8! 

1  ,  Q ) 

2  D( 

7  774 
1,0) 
3400 
3923 

i  ,  31 

■■,  ;o  • 

'♦  "':>  s 
1,0) 
I  6  ) 

8  133 
1  ,  Q ) 
700U 

3  3  78 
1,0) 

6  700 
1  9  J  8 
I,Q) 
3  200 
1  c^4? 

1  ,  )  > 

r,  -v  ) 

2  2  3  'i 
1,0 
5  0-1 

7  4  k 
It)) 
580- 
7  38  3 
1,0) 


*B 
2* 
4-FCT(  (  A-C)  ,  Z  ,0,T,EI  ,0)  ) 

*  B 
-2v 

+-PCT(  (  A-C)  ,Z,0,T,FI  ,C)  ) 
«8 

0-1* 

4-FCT  (  (A-C)  ,Z,D,T,EI,Q) ) 
*B 

D-l* 
4-FCT  (  (  A 
*B 
)-l* 
4-FCT  (  (A-C)  ,Z,U,TtEI,Q)) 

0-1* 

4-FCT  (  (  ft: 

*R 

0-1* 

+  FCT  (  (A- 

*B 

0-1* 

4-FCT  (  (  A- 

*B 

D-l* 

+-FCT  (  (  A- 

*  B 

D-l* 
4-FCT  (  (  A- 


.-C)  ,  Z,D,T,  F  I  ,0)  ) 


C)  ,Z,D,T, E I ,Q)  ) 


-F) , Z,D,T,EI,Q) ) 


.-C)  ,Z  ,0,T  ,  E  I  ,G)  ) 


C) ,Z,D,T,E l,Q) ) 


i-C)  ,Z,0,T,FI  ,0)  ) 


D-l* 

4  F  C  T  (  (  A 
*B 

0-1* 
4FCT( ( A 
*B 

D-l* 
4-FCT(  (  A 
■■-c\ 

0-1* 
4-FCT  (  {  A 

*B 
0- 1  * 

4-  F  c  r  (  (  A 
1  *  B 

g,  ;D-1* 
4-FCTl  (A 


-C),Z,D,T,EI,Q)) 


-C)  ,Z,0,T,FI  ,Q)  ) 


C),Z,D,T,EI,Q)) 


-C)  ,Z,3,T,EI , 3)  ) 


),Z,D,T,EI,Q)J 


-C)  , Z, 0,T, [  I , 0)  )  ) 
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c 
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c 
c 
c 
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REAL  FUNCTION  W I0N*8( X, EI , T , Z , A,Q  ) 

PURPOSE:  THIS  FUNCTION  CALCULATES  THAT  PART  OF  THE  ENERGY 
DISTRIBUTION  OF  AN  ELECTRON  HUE  !"U  LOSS  OF  ENERGY  BY 
IONIZATION  WHILE  PASSING  THROUGH  AN  ABSORBER. 

IMPLICIT  REAL*8( A-H,0-Z ) 

IL  ARE  THE  AVERAGE  IONIZATION  POTENTIALS  PER  ELECTRON  FOR 
LEAD  FOR  EACH  SHELL,  BEGINNING  WITH  THE  K  SHELL. 

NPOTL  IS  THE  NUMBER  OF  IONIZATION  POTENTIALS  USED  FOR  LEAD 

IN  THIS  CONTEXT,  L=LEAD,  G=GADOL I Nl UM,  S=TIN,  C=COPPER, 
A= ALUMINUM,  AND  B=3ERYLLlUM. 


C 
C 
C 

C 
C 
C 

C 
C 
C 

c 

C 


c 
c 
c 

c 
c 
c 


c 
c 
c 


)Ci  4  90 


, 


REAL*  8  IL('J  )/.  ;9  80U0,.  3 14700,. C(  320 
1  ,0<  DC  76,  .  '..   .  2/ 
■    R  E  AL*8  NL  (  6  )  / 2  .  ,  8  .  ,  1 8  .  ,  32 .  ,  18  .  ,  4.  / 

REAL*8  IG(5)/.wi>  24, .(  )785 , .00189 , .00024, .000026/ 

Re  AL*8  NG(5)/2.,3.,1S.,25.,9./ 

REAL* 8  IS (5  J/.C292,  .  v  44,  .bOC  72,  .000  )62,  .0C-C001/ 

REAL* 8  NS(  5)/2.  ,  -i*   ,  1°-.  ,1  3.  ,4./ 

REAL*8     IC(3)/.C0898,.(   3099, .000068/ 

REAL  !  6     NC(3 >/2. ,S. ,18./ 

RE AL*8  I A{  2  )  / .  (  • .  1  r>t   ,  . CO 1-088/ 

R!  AL*8  NA(  2)/2.  ,8./ 

Rl  \L*3  .18(2)/.   ■-■l*tf?,  .OOOC  14/. 

REAL* 3  NB(2)/2.,2./ 

i\  P  0  T  L  =  6 

NP0TG=5 

NP0TS=5 

N P OT  C= 3 

NP0TA=2 

Nj  P  0  y  B  =  ? 

SB  IS  THE  SUMMATION  OF  IONIZATION  POTENTIALS  FOR  THE 
EXPERIMENTAL  ATOM. 

SB=0. 
XMASE  IS  THE  REST  MASS  OF  AN  ELECTRON  IN  MEV. 

XMASE=. 511006 

f>FTA  IS  THE  NORMAL  V/C  USED  IN  RELATIVISTIC  CALCULATIONS, 

GAMMA  IS  THE  NORMAL  TERM  USED  IN  RELATIVISTIC  CALCULATIONS 

BETA=DSQRT( ( EI  *E I+2.*E I*XMASE ) / (E I*EI+2.*E I*XMASE 
1+XMASE*XMASE) ) 
GA 4MA=l./( L.-BETA**2)** (1./2. ) 

AT  THIS  POINT,  THE  PROGRAM  MUST  COMPUTE  THE  PROPER  CON- 
STANTS FOR  THE  EXPERIMENTAL  Z.   A  SERIES  OF  LOGIC  STATE- 
MENTS SELECT  THE  CORRECT  FORMULAE. 
SEE  \LUMINUM  AS  EXAMPLE. 

IE( Z.GT.4. 1 )G0T013 
DO  7  .  1  =  1,  JPGT8 

B=IB( I )*NB{ I )*QLQG{ 2.*EI/(  IB( I )*( l.-BETA**2) )  ) 
1       SB=S8+B 
82=18.64 

DELT=4.6C6*DL0G1  }{ 8ETA*GAMMA )-2.83 


DELT=4.6C6X+C  FCR  BERYLLIUM. 
GO  TO  111 
13  IF(Z.GT.13.1 )G0T02  9 

CCMPUTF  SB  I  II  AL  I 

TIAL  FROM  ONE  SHELL. 


SEE  STERNHEIMER. 


B  I S  THE  IONIZATION 


60 


c 
c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 
c 
c 
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c 

c 
c 
c 

c 
c 
c 
c 


DO  7  1  l=l,,MPOTA 

B=IAU)*NA(  I)*DLOG(2.*EI/(  IA(  I  )*(  l.-BETA**2))  ) 

71  SB=S3+B 

ESTABLISH  VALUES  FOR  32  AND  Al.   THESE  ARE  MATERIAL  DE- 
PENDENT CONSTANTS  USED  TO  COMPUTE  THE  AVERAGE  ENERGY 
LOSS  BY  IONIZATION  AND  ARE  SUPPLIED  BY  THE  STERNHLIMER 
PAPER.   THESE  ARE  B  AND  A,  RESPECT  I VELY f ' FROM  HIS  PAPER. 

82-16.77 
Al=.0740 

DELT  IS  A  PART  OF  THE  EXPRESSION  BY  STERNHEIMER  FOR  THE 
AVERAGE  ENERGY  LOSS  DUE  TO  IONIZATION.   NUMBERS  HEREIN 
REPRESENT  OTHER  STERNHEIMER  CONSTANTS  AS  FOLLOWS: 

SMALL  A   .J 906 

Sf^ALL  M   3.510 

XI        3. LUG 

-C        4.210 

DE  LT  =4. 606* DLOG 1 0 ( BETA*GAMMA J-4.21 +.0906* 
1  (3.-0L0GlC(RETA*GAMMA)  )**3.51 
GO  TO  111 
29  IFCZ.GT.29.  DGOT050 
DO  7  2  I=1,NP0TC 
B=IC(  I)*NC(I)*DLQG(2.*£I/{ICm*(l.-BETA**2n  ) 

72  SB=SB+R 
62=15.09 
Al=.0701 

DELT=4. 606*DLOG10( BE TA *GAMMA ) -4. 74  +  . 1 19C  * 
1 (3.-DLOG10(BETA*GAMMA) )**3.38 
GO    TO    111 
50     IE(Z.GT  .50.  DG0T064 
DO    7  3     I=1,NP0TS 
B=IS( I  )*NS( I )*DL0G(2.*EI/<  IS( I  )  * ( 1 . -BET A**2 ) )  ) 

73  SB=SB+B 
82=13.83 
Al=.r647 

DELT=4.6-:6*DL0G1C(  BE  TA*GAMM  A  )  -6  .  28  +  .  404*  (  3.- 
L0L0G10< BETA*GAMMA) ) 
1**2.52 

GO    TO    111 
64     IF(Z.GT.64.1)G0T082 

DO    74    I=1,NP0TG 

B=IG( I  )*NG<  I )*DLOG( 2.*EI/(  IG( I  )*( 1. -BET A** 2)  )  ) 

74  S8=SB+B 
B2=13.3 
Al=.0624 
DELT=4.606*DLOG10(BETA*GAMMA)-6.8+.418*(3.- 

1DLOG10(  dETAf-GAMMA)  ) 
1**2.10 
GC    TO    111 
82     IF(  Z.GT.82.DG0T0112 
DO    75     I=1,NP0TL 
B=IL( I )*NL( I )*DL0G(2.*EI/( IL( I )*( l.-BETA**2) ) ) 

75  S3=SR+B 
B2=12.81 
Al=.060  8 

DELT=4. 6C6*DL0G1Q( BETA*GAMMA ) -6.9 3+. 0652* 
1 ( 4.-DLOG10( BETA*GAMMA ) ) **3.41 
GO  TO  111 

AP  IS  THE  QUANTITY  SMALL  A,  FROM  BLUNCK  AND  WESTPHALt 
MULTIPLIED  BY  R,  THE  TARGET  THICKNESS  IN  CENTIMETERS. 

Ill  AR=0.154*Z*T/A/B ETA/BET A 

CCMPUTE  RS2,  THE  SMALL  B  SQUARED  FROM  THE  BLUNCK  AND 
LEISEGAfJG  PAPER. 

BS2=3.*SB/(Z*AR] 
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c 
c 

c 

c 

c 

c 
c 
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c 
c 
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c 
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c 
c 
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c 


c 
c 
c 


c 
c 
c 
c 
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c 
c 
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VAR1,  VAR2,  AND  VAR?  ARE  CONVENIENCE  STORAGE-  LOCATIONS 

USED  TC  STORE  COMPUTED  PORTIONS  OF  STERNHE I MER ■ S  AVERAGE 
ENERGY  LOSS. 

VAR1=A1*T/BETA**2 

VAR2=B2+.43  «-2.  *DLOG  (  BE  TA*GAMMA  ) 

VAR3=DLCG< El )-BETA**2 

QAVE  IS  TOT  FINAL  VALUE  FOR  THE  AVERAGE  ENERGY  LOSS. 

QAVE=VAR1* ( VAR2+VAR3-DELT) 
B1=DSQRT(8S2) 

X  IS  THE  PORTION  OF  SOME  TOTAL  ENERGY  LOSS  Q  DUE  TO 
I  ON  I Z AT  ION. 

X=Q-X 

XLAMB  IS  THE  LAMBDA  VARIABLE  FROM  THE  BLUNCK  AND  WESTPHAL 

P  A  P  :_  R  . 

XLAMB=< X-QAVE) /AR+DLOG ( EI/AR >-l • 1 16 
WION=0. 

TFF    VARIABLE    TEST     IS    USED    TO    RESTRICT    OUTPUT    TC    NUMBERS 
OF     LSEFUJL    MAGNITUDE. 

TEST=XLAM3**2/(B 1*81+3.24) 
IF( TE ST.G1  .  1  38  .  IGOTOl] 2 

COEFF    ARE    THE    COEFFICIENTS    FOR    EACH    OF    \Ht    FOUR    GAUSSIAN 
CURVES    BLUNCK    AND    WESTPHAL    USE    TO    FIT     IMF    LANDAU    DISTRI- 
BUTION.       AN    ADDITIONAL     CONSTANT     HAS    SEEN    MULTIPLIED     INTO 
EACH    FOR    CONVENIENCE     IN    RUNNING    THE    PROGRAM. 

r  n  ^FF 1  =  ( 0  .  1  74*  1  .  8  )  /  D  S  :>  <  1(81*81+1.8*1.8)     *  .  1 D  1 
COEFF2=(  :».     38*2.  -  )  /DS3RT(B1*R1  +  2.C  *2.J  )*. 107 
CUEFF3=(  J. CI  9* 3.     ) /OSQRT ( 81*81+3. 0*3.0 ) *. 1 D7 
C0EFF4=(  j.i      7*5.G)/DSQRT(B1*B1+5.U*5. -J)*.1D7 


GAUS  A -IE  THE  B- 


GAUSSIAN  CURVES  INCLUDING  COEFFICIENTS. 


GAUS1=C0EFF1*DEXP( - (XL AMB-C.O ) **2/ ( 81*81+1.3*1. 8) ) 
GAUS2=C0EFF2*Di  XP(  -(XL  '\MB-3.U  )**2/( B1*B 1+2.  i*2.  Z  )  ) 
GAUS3  =  C0EFF3*DEXP( -( XL  AM 3-6. 5 )** 2/ ( 81  *B 1  +  3.0*3. U)  ) 
GAUS4=C0EFF4*DEXP(  -(XL  XMB-ll.  ))**2/  (  Bl  *B1  +  5.<  *■>.  _  I  ) 

CCRINT  IS  \0  EMPERICALLV  DERIVED  CORRECTION  FACTOR  TO  I  IM- 
PENSATE  FOR  THE  -ACT  THAI  THE  D ISTR I3UT ION'  FUNCT ION  CAN- 
001  BE  INTEGRATED  EXACTLY  AS  MHE  RADIATION  LOSS  DISTRI- 
BUTION BECCMES  INFINITE  ^S  ITS  VARIABLE  (ENERGY  LOSS) 
APPROACHES  U\^% 


C0RINT=1.26*(1  .-.'. 
GAUS1=GAUS1*C0RINT 


a->  f*T**3.03) 


WION  IS  THE  SO:'  if  GAUSSIAN  CORVES  REPRESENTING  THE  ENERGY 
DISTRIBUTION  00c  TO  IONIZATION  LOSSES. 

W I ON=GAUS 1 +  0 AU  S 2 +0  AUS3 +GAU S  4 
1 12  RETURN 

TND 
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REAL  FUNCTION  WRAD*8( X ,E I tT ,Z f A , Q) 


PURPOSE:  THIS  FUNCT  1  IN 
DI  STR  IOUTICM  OF  AN  EL 
RACIATION  WHILE  PASS  I 


CALCULATES  THAT  PART  OF  THE  ENERGY 
ECTRON  DUE  TO  LOSS  OF  ENERGY  BY 
NG  THROUGH  AN  ABSORBER. 


IMPL  ICIT  REAL*8(  \-H,0-l) 

U  IS  AN  INTERMEDIATE  STORAGE  LOCATION  USED  IN  COMPUTING 
ALPHA. 


U=183.C/Z#*< 1. /3.0 ) 

COMPUTE  VALUE  OF  ALPHA' 
NUMBEI  AND  ATOM  IC  -.'E  I 
THE  TARGET  THICKNESS 

ALPHR=u. j014*T*Z*Z/ 

CCMPUTE  VALUE  FOR  B,  A 

B=1./DGAMMA< ALPHR+1 

COMPUTE  WRAD,  THE  MUMBE 
RAD  I  AT  ION  ENERGY  LO.SS 
SLVL'  TOTAL  ENtP.GY  L  jS 


R  FOR  TARGET  THICKNESS,  ATOMIC 

GHT  OF  EXPERIMENTAL  ABSCR3EP.   R  IS 

EXPRESSED  IN  CENT  I  "~  T  h  -<S  . 

A*(4.O/3.C*DLOG(U)+1.0/<?.L) 

NORMALIZING  CONSTANT. 

6  ) 

R  OF  COUNTS  EXPECTED  AS  A  RESULT  OF 
ES.   X  REPRESENTS  THAT  PORTION  OF 
S  WHICH  IS  NOT  LOST  BY  IONIZATION. 


W R  A  D= ( B *  A L  P HR * ( X  -  *  A  LPHR ) / ( E I ** ALPHR ) ) /X 

RETURN 
END 


C 
C 

c 

c 
c 
c 
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REAL  FUNCTION  BE AM*8( E I , HW, X ) 

PURPOSE:  THIS  FUNCTION  REPLICATES  THE  HSTRIBUTIGN  OF 
BEAM  ENERGY  AS  A  S I  "•)  ^  L r:  GAUSSIAN  CURVE  BASED  ON  ! 
EXPERIMENTALLY  OBSERVED  HALFfolDTH  AND  IS  LOCATED 
RELATIVE  TO  THE  BEAM  ENERGY  BY  THE  COMPUTED  AVt;^AGE 
RECOIL  LOSS. 

IMPL ICI r  RE AL*8( A-H,0-Z) 

Y=2. 

ALFA=4.*DL0G<Y  )/( HW*HWJ 

BEAM=DI  XP{-{  (X-E  I  )  **2)  -ALFA) 

RETURN 

END 


REAL  FUNCTION  FCT*8 ( XT Z , A , T t E I , 0) 
C 

C  PURPOSE:  THIS  FUNCTION  "ULT I  PL  I  CAT  I VF LY  JOINS  TOGETHER  THE 
C    PREDICTION  OF  ENERGY  LOSS  DUE  TO  IONIZATION  AND  THE  LOSS 
C    DUE  TO  RAD  I  AT  ION. 
C 

IMPL ICIT  RFAL*8( A-HtO-Z) 

FCT=WRAD(X,EI  ,  T,Z,  A,Q     )  •■  .   I      KX,EI»TfZtAtG     ) 
;  TURN 

END 
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SUB  ROUT  I 'ME  GRAPH  (  THX  ,  THY  t  EXX  ,  EXY  ,  T  1  ,  T2  ,  T3  »  T4  t  HW  ,  E  I  ♦  T  , 
1C0RR,CTSNQR,QP) 

PURPOSE:  THIS  SUBROUTINE  PLUTS  THE  RESULTS  OF  THEORY  CAL- 
CULATIONS AND  EXPERIMENTAL  RESULTS. 

IMPLICI T  REAL* 8 ( A-H,Q-Z) 
REAL*d  T  ITLEX(  1 ) 
C  A  T  A  T  I T  L  E  X  / '  M  if  V      «  / 
REAL*8  TITLI  Y( 3) 

DATA  TITLEY/'NORMALIZEO  COUNTS        •/ 
REAL* 8  TITLEU(2J 

CATA  TITLEG/'QPPEDAHL  BOX  0   • / 
D I  ME"  \S  I  ON  T 1 < 4 ) , T2 ( 4 ) , T3 <  4 ) , T4 ( 4 ) , THX ( 3C 
1EXX(20  )  ,EXY(?-  J)  ,  CTSUP(2UC  )  fCTSDN(20  ,) 


)  ,THY(  3CC  )  , 


DETERMINE    SCALING    FACTORS     AND    CHANGE    THEORETICAL    CURVE 
TO    PRO'ER    LOCATION. 

X  =  f<. 

Y  =  5. 

SIXIN=( E l-OP)*2. 

ISIXIN=SIXIN+.  5 

SIXIN2=ISIXIN 

SI XIN=SIXlN2/2. 

IF(HW.EQ.O)GOT  )90 

XNCR2=i5*HW 

INCR2=XNCR24 .5 

IF(  INCR2.LT. 1 )  I  NCR 2  =  1 

XINCR2= INCR2 

XINCR=XINCR2/2. 

SP  =  SIXlN-6,-  XI  NCR 

[F( SP.LT.O )G0T09C 

GOT 091 

90  SP  =  C 

INCR  =  (  EI-1.5*T  1/6.  +  .5 
X1NCR=INCR 

91  N  =  '., 

2C  1=1+1 

IE  (E XX (  I  ).GT.l.  )  G0T021 

GOTO 22 

N=N  +  J 

GOTO  20 

IFCN.EO.yjGOTGlO 

DO    23    J-lfN 

EXX(J)=(EXX{J)-SP) /XINCR+CORR/XINCR 

IF{ J. E Q.N) GOTO 23       - 
23    CONTINUE 

I  =  J 

•v  =  I 

1  =  1  +  1 

1F(THX( I ) .GT 

GOT022 

M=M+1 

GOTO  30 

DO    33    J=1,M 

THX (J) =(THX( J)-SP) / XI  NCR 

IF( J.EQ.MIG0T03  3 

IF(THY( J+1).GT.THY< J) ) ThYMAX=THY( J+ 1 ) 
33    CONTINUE 

IF(N.ECUG)G0TO101 

SCALE=4. 

IF(T.GT .4.  3)     SCALE=3. 

DO    4  )    J  =  1,N 

LXY( J) =EXY( J )/CTSNOR*SCALE 
A.:    CO  IT!  YJE 
i     1     DO    5       J  =  1  »   ' 

SCALE=4. 

IF{T.GT.4.8)     SCALE=3. 


21 
22 


100 
30 

31 

32 


1. )     GOTO    31 
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c 
c 
c 


c 
c 

c 


THY( J)=THY< J >/THYMAX*SCALE 

50    CCMTINUE 

SORT     EXX    AND    ELIMINATE    VALUES    OFF    GRAPH. 

IF(N.EQ.O) GOTO  102 

K=N-1 

DC    60     I =  1  ,  K 

IP  1=1  +  1 

DG    60     J=IPltM 

IF(E>X(I).LL".EXX(J))    GOTO    6u 

TEMP=EXX( I ) 

TEM=EXY( I ) 

EXX( l)=EXX( Jl 

EXY(  1  )  =  EXY( J ) 

EXX(J)=TEMP 

EXY( J)=TEM 

60  CONTINUE 
J=J 

KK  =  0 
K  =  N 

DO    61     1=1, K 
IF(KK.F0.1 >GG1    161 
IF(EXX(I).LT.O  )GUT06l 
J  =  J  +  3 
N  =  J 

EXX( J)=EXX( 1 ) 
EYY< J)=EXY( I ) 
I  F ( GXX( J ) .GT.X )GOT062 
GOTO    61 
62    N=J-1 
KK=1 

61  CONTINUE 
102  KK=0 

J  -  ... 

i<  =M 

DC  71  I =1 »K 

IF(KK.EQ.  IJGOT071 

IF(TnX(  I  ).LT.0 )GOT071 

J  =  J+1 

M  =  J 

THX( J)  =  THX(  1  ) 

ThY( J)=THY( I ) 

IF( THX(  I  ) .01 '. X )GOT  372 

G0TO71 
72  M=J-1 

KK  =  1 
71  CONTINUE 

DETERMINE  ERROR  BAP  MAGNITUDE, 

IF(N.EQ.O) GGT0103 
DO  7  0  I  =1,N 

CTSDM( I )=bXY( I |-DSQRT(EXY( I )*200. )/2.0. 
CTSUP( I ) =EXY( I >+USQRT(EXY(I >*2C0. )/2  -. 
700  CONTINUE 

INITIALIZE  PLOT  AND  WRITE  IDENTIFICATION. 

1  03  CALL  PI  ITS 

CALL  SYM30L(  t  ,  .28 ,T ITLEO.G , 16 ) 
CALL  PLOTi : i 2.t-3) 

ORAw  OUTLINE  OF  GRAPH. 

Z=0 

CALL  PL0T(Z,Y,2J 

CALL  PL0T(X,Yt2) 

CALL  PL0T(X,Z,2) 

CALL  PLOT (Z ,Z, 2) 

DRAW  OUTLINE  OF  TITLE  BCX. 
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7  U 


Bl 
82 
B3 
B4 
B5 
8  6 
CA 
CA 
CA 
CA 
CA 
CA 
B5 
CA 
B5 
Bo 
EX 
CA 
R5 
ob 
CA 
3b 
86 
CA 
IF 
86 
3  5 
B7 
CA 
CA 
B6 
B5 

r>n 

CA 

Rn 
r  a 


.4 
4.4 

3.  54 
4.6 

:  B 1  +  . 

=  B3+. 


PL 
PL 
PL 
PL 
PL 
SY 
5  +  . 
NJ 
3  5+1 
B6  +  . 
=  2. 
L  NU 
85  +  1 
86- . 
L  W 
B5-2 
B6+. 
L  SY 
N,  fc  Q 
B6+. 
'1  5  +  * 
B5  +  . 
L  PL 
L  PL 
3  5-. 
8  5-. 
36+. 
L  SY 
R  6  +  . 
L    SY 


T 1 ,  0  •  G ,  32  ) 


1 
1 

0  T  ( B 1  ,  B  3  ,  3 ) 
]T(32,B3, 2) 
0T(B2»B4 ,2) 
OT(Bl,B4,2) 
0T(B1, B3,2) 
MBOL( B5 , B6, .  14  T 
24 

M8ER(B5  ,  B6,  »  14, T  ,0.8 , 3  ) 
.  32 
1 

MBER(B5,B6,.07t£XP,U.C,-l) 
.20 

i 

MBER(B5t Bbt.  14, EI  ,0. 3 , 2 ) 

.76 

24 

M80L(B5 , B6,. 14, 

•  o) GOTO  7    I 

J  7     .   / 

P2      '      • 

58 

0T-    (B5,36,3) 

0  T     (87,36,2) 

0  7 

82 


M.80L(B5»  8  6,  .  14, 

?4 

M.8QL  (  B5»  86,  •  14, 


C 

C    PLOT    TIC    "MRKS. 

C 

G  =  .  1 

H  =  Y—     1 

15  7U    C'flLL*PL0T(G,Y,3J 

CALL    PLOT(G,H,2) 

G=G+. 1 

IF     (G.LE.X)GO    TO    151  j 

G  =  X-. 1 

H=.l 
158C    CALL    PLOT(G,H,3) 

CALL    PLOT(G,0.,2) 

G  =  8—    1 

IF(G*GT.G. )G0    TO    15  30 

G  =  l. 

H=.2 
167.,    CALL     PLOTCGtH,  3) 

CALL     PLOT(G,, , 2) 

G  =  G+ 1 

fpfG.LT.X)     GO    TO    167:. 

H  =  Y~     2 

16  88    CALL* PLOT (G,H, 3) 

CALL     PLJT(G,Y, 2) 

C = C  —  1 

IF(G.GT.( . )     GO    TO    168 ) 

G=.l 

u—  y_     i 

1571     CALL* PLOT (G,H, 3) 

CALL     PLOT (_.,H, 2) 
H=H-.l 

IF(H.GT.  ■•.  )     GO    TO    1  57  1 
G  =  X-. 1 
H=.  1 
1581     CALL     PI  0T(  3  ,  H  ,  3  ) 
CALL     PL0T(X,H,2) 


T3,0.<   ,  32) 
T4 , G •  I   ,32) 
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H  =  H  +  .  1 

IF(H.LT.Y)  GU  TO  1581 

G=X-.2 

H=Y-1. 

16  71  CUL  PL0T(oiH,3) 
CALL  PLUT(X,H,2) 

H  —  H—  1 

IFtH.GT. ))  GO  TO  1671 
G=.2 

H=l. 
1681  CALL  PL0T<G,H,3) 
CALL  PLOT(  ,H, ?  ) 
H=H+1 . 
I F ( H . LT • Y I  G i J  TO  1681 

LAaEL  AXIS. 


C 

r 

c 


CALL  SY  -1BOH  3  .r>,~.  6  ,.  14  .TITLEX,  [  ,4) 
CALL  SYMBOL (-.2 5, i . 56 , . 14  f TITLEY f 9 »., 24) 

NUMBER  ENERGY  AXIS. 

10  50    G=-.2  8 

H=-.32 

FLO=SP 

1C6C     CALL    NUMBER (G,H, .14tFL0»0,2) 

G=G+1, 

IF(G.GT.X)  GO  T  ]  17 
FLO=FLO+XI NCR 
Gj  TO  lu60 
C 

C  PLOT  EXPERIMENTAL  DATA. 
C 
1,7   IF(N.E  ).0)G0T01C4 

r M  L  L I NE ( t XX,EXYf N»2, -6) 

r 

C  PLCT  ERROR  FLAGS. 
C 

DO  1.1   1  =  1  fN 

ES1=FXX(I)-.  5 

ES2=EXX( I )+.  b 

CALL  PL;1T(L  Sl,C  TSUP(  I)  t3) 

CALL    PLUT(ES2i CTSUP( I )  , 2  ) 

CALL     PLOT (bXX ( I  )  ,CT  SUP(  I  )  ,3  ) 

CALL     PLOT(EXX(  I  )  ,C  I  SOM<  I)  ,2) 

CALL     PL'JT(  F  S]  ,  CTSD'N(  1  )  ,  3  ) 

CALL    PL  OT(!  S2,CTSDN( I)  ,?  ) 
i     10    CONTINUE 
C 

C    PLCT    THEORY    CUTU'l  . 
C 

104    CALL    LINElTHXf THYfMf2tl) 
C 

C    F  INAL  IZE     PLOT. 
C 

CALL    PLOT(   i.,8.,-3) 

CALL    PLOTE 

RETURN 

END 
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